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Abstract. We study enhancement of diffusive mixing on a compact Riemannian manifold 
by a fast incompressible flow. Our main result is a sharp description of the class of flows that 
make the deviation of the solution from its average arbitrarily small in an arbitrarily short 
time, provided that the flow amplitude is large enough. The necessary and sufficient condition 
on such flows is expressed naturally in terms of the spectral properties of the dynamical 
system associated with the flow. In particular, we find that weakly mixing flows always 
enhance dissipation in this sense. The proofs are based on a general criterion for the decay 
of the semigroup generated by an operator of the form F + iAL with a negative unbounded 
self-adjoint operator F, a self-adjoint operator L, and parameter A ^ I. In particular, 
they employ the RAGE theorem describing evolution of a quantum state belonging to the 
continuous spectral subspace of the hamiltonian (related to a classical theorem of Wiener on 
Fourier transforms of measures). Applications to quenching in reaction-diffusion equations 
are also considered. 



Let M be a smooth compact li-dimensional Riemannian manifold. The main objective 
of this paper is the study of the effect of a strong incompressible flow on diffusion on M. 
Namely, we consider solutions of the passive scalar equation 



Here A is the Laplace-Beltrami operator on M, m is a divergence free vector field, V is the 
covariant derivative, and A G M is a parameter regulating the strength of the flow. We are 
interested in the behavior of solutions of (jl.lj) for A 3> 1 at a fixed time r > 0. 

It is well known that as time tends to infinity, the solution t) will tend to its average. 



with \M\ being the volume of M. We would like to understand how the speed of convergence 
to the average depends on the properties of the flow and determine which flows are efficient 
in enhancing the relaxation process. 

The question of the influence of advection on diffusion is very natural and physically 
relevant, and the subject has a long history. The passive scalar model is one of the most 
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studied PDEs in both mathematical and physical literature. One important direction of 
research focused on homogenization, where in a long time-large propagation distance limit 
the solution of a passive advection-diffusion equation converges to a solution of an effective 
diffusion equation. Then one is interested in the dependence of the diffusion coefficient on 
the strength of the fluid flow. We refer to [21] for more details and references. The main 
difference with the present work is that here we are interested in the flow effect in a finite 
time without the long time limit. 

On the other hand, the Freidlin- Wentzell theory [121 El UHl HHI studies in and, for 
a class of Hamiltonian flows, proves the convergence of solutions as A ^ oo to solutions of 
an effective diffusion equation on the Reeb graph of the hamiltonian. The graph, essentially, 
is obtained by identifying all points on any streamline. The conditions on the flows for 
which the procedure can be carried out are given in terms of certain non- degeneracy and 
growth assumptions on the stream function. The Freidlin- Wentzell method does not apply, 
in particular, to ergodic flows or in odd dimensions. 

Perhaps the closest to our setting is the work of Kifer and more recently a result of Beresty- 
cki, Hamel and Nadirashvili. Kifer's work (see [211 122 1211 1211 where further references can 
be found) employs probabilistic methods and is focused, in particular, on the estimates of 
the principal eigenvalue (and, in some special situations, other eigenvalues) of the operator 
— eA + M ■ V when e is small, mainly in the case of the Dirichlet boundary conditions. In 
particular, the asymptotic behavior of the principal eigenvalue Aq and the corresponding pos- 
itive eigenfunction 0q for small e has been described in the case where the operator u ■ V 
has a discrete spectrum and sufficiently smooth eigenfunctions. It is well known that the 
principal eigenvalue determines the asymptotic rate of decay of the solutions of the initial 
value problem, namely 

limt-Mog||0^(x,t)|U2 = -A^ (1.2) 

(see e.g. [221)- In a related recent work [21, Berestycki, Hamel and Nadirashvili utilize PDE 
methods to prove a sharp result on the behavior of the principal eigenvalue of the operator 
—A + Au ■ V defined on a bounded domain f2 C M'^ with the Dirichlet boundary conditions. 
The main conclusion is that A^ stays bounded as A — >■ oo if and only if u has a first integral 
w in HI{VL) (that is, u ■ Vw = 0). An elegant variational principle determining the limit of 
Aa as A — > oo is also proved. In addition, [21 provides a direct link between the behavior of 
the principal eigenvalue and the dynamics which is more robust than ()1.2|) : it is shown that 
||0"^(-, 1)11^2(^2) can be made arbitrarily small for any initial datum by increasing A if and only 
if A^ ^ cxD as y4 ^ oo (and, therefore, if and only if the flow u does not have a first integral 
in Hq{VL)). We should mention that there are many earlier works providing variational 
characterization of the principal eigenvalues, and refer to [21 123 for more references. 

Many of the studies mentioned above also apply in the case of a compact manifold with- 
out boundary or Neumann boundary conditions, which is the primary focus of this paper. 
However, in this case the principal eigenvalue is simply zero and corresponds to the constant 
eigenfunction. Instead one is interested in the speed of convergence of the solution to its 
average, the relaxation speed. A recent work of Franke provides estimates on the heat 
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kernels corresponding to the incompressible drift and diffusion on manifolds, but these esti- 
mates lead to upper bounds on ||</''^(1) — 0|| which essentially do not improve as A ^ oo. 
One way to study the convergence speed is to estimate the spectral gap - the difference 
between the principal eigenvalue and the real part of the next eigenvalue. To the best of 
our knowledge, there is very little known about such estimates in the context of (jl.lj) : see 
[221 P- 251 for a discussion. Neither probabilistic methods nor PDE methods of j2] seem to 
apply in this situation, in particular because the eigenfunction corresponding to the eigen- 
value(s) with the second smallest real part is no longer positive and the eigenvalue itself does 
not need to be real. Moreover, even if the spectral gap estimate were available, generally it 
only yields a limited asymptotic in time dynamical information of type (jl.2p . and how fast 
the long time limit is achieved may depend on A. Part of our motivation for studying the 
advection-enhanced diffusion comes from the applications to quenching in reaction-diffusion 
equations (see e.g. jH [T21 123 EH EH] ) , which we discuss in Section [7| For these applications, 
one needs estimates on the A-dependent L°° norm decay at a fixed positive time, the type 
of information the bound like ()1.2|) does not provide. We are aware of only one case where 
enhanced relaxation estimates of this kind are available. It is the recent work of Fannjiang, 
Nonnemacher and Wolowski [101 E], where such estimates are provided in the discrete setting 
(see also [22] for some related earlier references). In these papers a unitary evolution step (a 
certain measure preserving map on the torus) alternates with a dissipation step, which, for 
example, acts simply by multiplying the Fourier coefficients by damping factors. The absence 
of sufficiently regular eigenfunctions appears as a key for the lack of enhanced relaxation in 
this particular class of dynamical systems. In [IHIIIII, the authors also provide finer estimates 
of the dissipation time for particular classes of toral automorphisms (that is, they estimate 
how many steps are needed to reduce the norm of the solution by a factor of two if the 
diffusion strength is e). 

Our main goal in this paper is to provide a sharp characterization of incompressible flows 
that are relaxation enhancing, in a quite general setup. We work directly with dynamical 
estimates, and do not discuss the spectral gap. The following natural definition will be used 
in this paper as a measure of the flow efficiency in improving the solution relaxation. 

Definition 1.1. Let M be a smooth compact Riemannian manifold. The incompressible flow 
u on M is called relaxation enhancing if for every r > and 5 > 0, there exist yl(r, 5) such 
that for any A > A{t, 5) and any G L'^{M) with ||0o||l2(a/) = ^ we have 

U%,t)-^lhm)<S, (1.3) 

where (j)"^{x,t) is the solution of (jl.lj) and the average of (pQ. 

Remarks. 1. In Theorem 15.51 we show that the choice of the norm in the definition is not 
essential and can be replaced by any L^-norm with 1 < p < oo. 

2. It follows from the proofs of our main results that the relaxation enhancing class is not 
changed even when we allow the flow strength that ensures (jl.3|) to depend on 0o, that is, if 
we require p.3j) to hold for all 0o G L^(M) with ||0o||L2(Af) = 1 and all A > yl(r, 5, 0o). 

Our first result is as follows. 
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Theorem 1.2. Let M be a smooth compact Riemannian manifold. A Lipschitz continuous 
incompressible flow u G Lip(M) is relaxation enhancing if and only if the operator u ■ V has 
no eigenfunctions in H^{M), other than the constant function. 

Any incompressible flow u G Lip(M) generates a unitary evolution group W on L'^{M), 
defined by U^f{x) = Here is a measure preserving transformation associated 

with the flow, defined by = M(<l>t(x)), ^o{x) = x. Recall that a flow u is called weakly 

mixing if the corresponding operator U has only continuous spectrum. The weakly mixing 
flows are ergodic, but not necessarily mixing (see e.g. 0). There exist fairly explicit examples 
of weakly mixing flows P HSl CH I2H1 ISHl IBS] , some of which we will discuss in Section El A 
direct consequence of Theorem 11.21 is the following Corollary. 

Corollary 1.3. Any weakly mixing incompressible flow u G Lip(M) is relaxation enhancing. 

Theorem II. 2t as we will see in Section |3J in its turn follows from quite general abstract 
criterion, which we are now going to describe. Let F be a self-adjoint, positive, unbounded 
operator with a discrete spectrum on a separable Hilbert space H. Let < Ai < A2 < . . . 
be the eigenvalues of F, and cj the corresponding orthonormal eigenvectors forming a basis 
in H. The (homogenous) Sobolev space if'"(F) associated with F is formed by all vectors 
ip = CjCj such that 

Note that H'^{T) is the domain -D(F) of F. Let L be a self-adjoint operator such that, for 
any ip G -f^^(F) and t > we have 

< C\mHHT) and lle^^VllHHD < 5(t) ||^|ki(r) (1-4) 

with both the constant C and the function B(t) < 00 independent of ip and B{t) G Lf^ci^, 00). 
Here e*^* is the unitary evolution group generated by the self-adjoint operator L. One might 
ask whether one of the two conditions in (jl.4j) does not imply the other. We show at the end 
of Section |21 by means of an example, that this is not the case in general. 
Consider a solution (p'^it) of the Bochner differential equation 

J^^\t) = ^AL<P^it)-T^\t), 0^(0) = 00. (1.5) 

Theorem 1.4. Let T be a self-adjoint, positive, unbounded operator with a discrete spectrum 
and let a self-adjoint operator L satisfy conditions p.4j) . Then the following two statements 
are equivalent: 

• For any t,6 > there exists A{t, 6) such that for any A > A{t, 6) and any (pQ E H 
with ||0o||// = 1; the solution (p^(t) of the equation (|1.5p satisfies ||0'^(t)|||^ < 6. 

• The operator L has no eigenvectors lying in H^{r). 

Remark. Here L corresponds to m ■ V (or, to be precise, a self-adjoint operator generating 
the unitary evolution group [/* which is equal to iu-V on H^{M)), and F to — A in Theorem 
O with H C L^{M) the subspace of mean zero functions. 
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Theorem 11.41 provides a sharp answer to the general question of when a combination of 
fast unitary evolution and dissipation produces a significantly stronger dissipative effect than 
dissipation alone. It can be useful in any model describing a physical situation which involves 
fast unitary dynamics with dissipation (or, equivalently, unitary dynamics with weak dissi- 
pation). We prove Theorem 1 1.41 in Sectional The proof uses ideas from quantum dynamics, 
in particularly the RAGE theorem (see e.g., [B|) describing evolution of a quantum state 
belonging to the continuous spectral subspace of a self-adjoint operator. 

A natural concern is if the existence of rough eigenvectors of L is consistent with the con- 
dition ()1.4|) which says that the dynamics corresponding to L preserves H^{T). In Section 0] 
we answer this question in the affirmative by providing examples where rough eigenfunctions 
exist yet (jl.4|) holds. One of them involves a discrete version of the celebrated Wigner-von 
Neumann construction of an imbedded eigenvalue of a Schrodinger operator Moreover, 
in Section ini we describe an example of a smooth flow on the two dimensional torus with 
discrete spectrum and rough (not if^(T^)) eigenfunctions - this example essentially goes back 
to Kolmogorov J8^. Thus, the result of Theorem 11.41 is precise. 

In Section [71 we discuss the application of Theorem 11.21 to quenching for react ion- diffusion 
equations on compact manifolds and domains. This corresponds to adding a non-negative 
reaction term /(T) on the right hand side of (jl.ll) . with /(O) = /(I) = 0. Then the long- 
term dynamics can lead to two outcomes: ^ 1 at every point (complete combustion), or 
(p^ c < 1 (quenching). The latter case is only possible if / is of the ignition type, that 
is, there exists 6q such that /(T) = for T < 6q, and c < 6*0. The question is then how the 
presence of strong fluid flow may aid the quenching process. We note that quenching/front 
propagation in infinite domains is also of considerable interest. Theorem II. 2l has applications 
in that setting as well, but they will be considered elsewhere. 

2. Preliminaries 

In this section we collect some elementary facts and estimates for the equation ()1.5p . 
Henceforth we are going to denote the standard norm in the Hilbert space if by || ■ ||, the 
inner product in H by (■, ■), the Sobolev spaces W^iV) simply by if™ and norms in these 
Sobolev spaces by || ■ H^. We have the following existence and uniqueness theorem. 

Theorem 2.1. Assume that for any ip G H^, we have 

IIL^II <C||^||i. (2.1) 
Then for any T > 0, there exists a unique solution (f){t) of the equation 

0'(t) = (zL-r)0(t), m = ^oeH\ 

This solution satisfies 

0(t) G L\[0,T],H')nC{[0,T],H'), <P\t) G ^^([0, T], ii). (2.2) 

Remarks. 1. The proof of Theorem l2.1l is standard, and can proceed by constructing a weak 
solution using Galerkin approximations and then establishing uniqueness and regularity. We 
refer, for example, to Evans [S] where the construction is carried out for parabolic PDEs but, 
given the assumption ()2.1|) . can be applied verbatim in the general case. 
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2. The existence theorem is also vahd for initial data 0o ^ but the solution has rougher 
properties at intervals containing t = 0, namely 

e L\%T],H^)nC{\Q,TlH), <P'(t) G L\[Q,T],H-'). (2.3) 

The existence of a rougher solution can also be derived from the general semigroup theory, by 
checking that iL — T satisfies the conditions of the Hille-Yosida theorem and thus generates 
a strongly continuous contraction semigroup in H (see, e.g. [Zj). 

Next we establish a few properties that are more specific to our particular problem. It will 
be more convenient for us, in terms of notation, to work with an equivalent reformulation of 
()1.5|) . by setting e = and rescaling time by the factor e~^, thus arriving at the equation 

m{t) = {tL - er)0^(t), 0^(0) = 00. (2.4) 

Lemma 2.2. Assume ()2.1|) . then for any initial data 0o £ H, \\4>o\\ = 1, the solution 4>''{t) 
of ()2.4|) satisfies 








Proof Recall that if G H\T), then T(f) G H-\T) and (10,0) = ||0||?. The regularity 
conditions ()2.2|l - ()2.3j) and the fact that L is self-adjoint allow us to compute 

|il01P = (0^0^) + (0.^0^) = -2e||01|^ (2.6) 
Integrating in time and taking into account the normalization of 0o, we obtain ()2.5j] . □ 

An immediate consequence of ()2.6|) is the following result, that we state here as a separate 
lemma for convenience. 

Lemma 2.3. Suppose that for all times t G {a,b) we have \\4>'^{t)\\l > A^||0^(t)|p. Then the 
following decay estimate holds: 

||0^(&)f <e-2^^(^-'^)||0^(a)f. 

Next we need an estimate on the growth of the difference between solutions corresponding 
to e > and e = in the if-norm. 

Lemma 2.4. Assume, in addition to ()2.1|) . that for any ■0 ^ o-nd t > we have 

||e^^Vl|i<5WII^I|i (2.7) 
for some B{t) < oo such that B{t) G Lf^^[0,oo). Let 0'^(t), 0'^(t) be solutions of 

(0°)'(t) = ^L0°(t), (0^)'(t) = i^L - er)0^(t), 
satisfying 0°(O) = 0*^(0) = 0o G H^. Then we have 

Remark. Note that 0°(t) = e*^*0o by definition. Assumption ()2.7p says that this unitary 
evolution is bounded in the H^{r) norm. 
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Proof. The regularity guaranteed by conditions ()2.H) . ()2.7|) and Theorem 12.11 allows us to 
multiply the equation 

(0^ - 00)' = iL{(f)' - 0°) - er0^ 

by 0*^ — 0". We obtain 

|||0^ - 0°f < 2e(l|0l|il|0°||i - II01I?) < lel|0°||?, 

which is the first inequality in (|2.8p . The second inequality follows simply from the assump- 
tion (EZZI). □ 

The following corollary is immediate. 

Corollary 2.5. Assume that fl2.1|) and ()2.7|) are satisfied, and the initial data 0o € H^. Then 
the solutions (t)''{t) and (jPit) defined in LemmaW^ satisfy 



\W{t)-<p\t)r<\e\\U\f^B\t)dt 

for any time t < t. 

Finally, we observe that conditions 1)2.11) and ()2.7p are independent. Taking L = T shows 
that ()2.7|) does not imply ()2.H) . because in this case the evolution e*'^* is unitary on but 
the domain of L is C H^. On the other hand, 1)2.11) does not imply ()2.7j) . as is the case in 
the following example. Let H = -L^(0, 1), define the operator F by F/(x) = e" /(n)e^™^ 
for all / G -ff such that e"^/(n) G ^^(Z), and take Lf{x) = xf{x). Then L is bounded on if 
and so ()2.1|) holds automatically, but 

so that e^'^^-^e^'^*"'^ = e^'^*("+^)'^. It follows that e^'"^^ is not bounded on (and neither is 
e*-^* for any t 7^ 0). 

3. The Abstract Criterion 

One direction in the proof of Theorem 11.41 is much easier. We start by proving this easy 
direction: that existence of H^iT) eigenvectors of L ensures existence of r, 5 > and 0o 
with 1 1 00 1 1 = 1 such that ||0'^(t)|| > 5 for all A - that is, if such eigenvectors exist, then the 
operator L is not relaxation enhancing. 

Proof of the first part of Theorem \1.4\ Assume that the initial data 0o G for ()1.5|) is an 
eigenvector of L corresponding to an eigenvalue E, normalized so that ||0o|| = 1. Take the 
inner product of p.5|) with 0o. We arrive at 

4(0-^(t),0o) =^AE(0-^(t),0o) - (F0^(t),0o). 



dt 

This and the assumption 0o G lead to 

|(e-^^*(0^(t),0o)) 



<l{u^mi + \\Mi) 
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Note that the value of the expression being differentiated on the left hand side is equal to 
one at t = 0. By Lemma 12.21 (with a simple time rescaling) we have \\(j)^{t)\\1dt < 1/2. 
Therefore, for t < r = (2||0o||?)"^ we have |(0^(t),0o)| > 1/2. Thus, ||0^(r)|| > 1/2, 



Note also that we have proved that in the presence of an eigenvector of L, enhanced 
relaxation does not happen for some 0o even if we allow A{t, 6) to be 0o-dependent as well. 
This explains Remark 2 after Definition 11.11 

The proof of the converse is more subtle, and will require some preparation. We switch to 
the equivalent formulation (j2.4j) . We need to show that if L has no eigenvectors, then 
for all T,6 > there exists eo(T, 5) > such that if e < eo, then ||0''(T/e)|| < 6 whenever 
1 1 00 1 1 = 1- The main idea of the proof can be naively described as follows. If the operator 
L has purely continuous spectrum or its eigenfunctions are rough then the if^-norm of the 
free evolution (with e = 0) is large most of the time. However, the mechanism of this effect 
is quite different for the continuous and point spectra. On the other hand, we will show that 
for small e the full evolution is close to the free evolution for a sufficiently long time. This 
clearly leads to dissipation enhancement. 

The first ingredient that we need to recall is the so-called RAGE theorem. 

Theorem 3.1 (RAGE). Let L be a self-adjoint operator in a Hilbert space H. Let he the 
spectral projection on its continuous spectral suhspace. Let C he any compact operator. Then 
for any 0o ^ H, we have 



Clearly, the result can be equivalently stated for a unitary operator [/, replacing e*^* with 
f/*. The proof of the RAGE theorem can be found, for example, in jH]. 

A direct consequence of the RAGE theorem is the following lemma. Recall that we denote 
by < Ai < A2 < . . . the eigenvalues of the operator V and by ei, 62, . . . the corresponding 
orthonormal eigenvectors. Let us also denote by Pn the orthogonal projection on the subspace 
spanned by the first N eigenvectors ei, . . . , ca? and hj S = {(p E H : ||</>|| = 1} the unit sphere 
in H. The following lemma shows that if the initial data lies in the continuous spectrum of 
L then the L-evolution will spend most of time in the higher modes of F. 

Lemma 3.2. Let K G S he a compact set. For any N,(t > 0, there exists Tc{N, a, K) such 
that for all T > Tc{N, a, K) and any (p & K, we have 



uniformly in A. 



□ 



T 








T 




(3.1) 







Remark. The key observation of Lemma 13.21 is that the time Tc{N,a, K) is uniform for all 
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Proof. Since Pjy is compact, we see that for any vector (j) E S, there exists a time Tc{N, a, (p) 
that depends on the function such that (j3.H) holds for T > Tc{N, a, (p) - this is assured by 
Theorem IH.ll To prove the uniformity in 0, note that the function 

T 



is uniformly continuous on S for all T (with constants independent of T): 

T 

|/(T,0) - /(T,^)| < |||P^e^^*Pe0|| - ||PAre^^*Pc^||| (||Pive^"Pc0|| + ||P^e^^*Pe^||) cit 



T 

<(II</'IH-Il^ll)^ / ||P^e'"*Pc(0-^)IMt< 2110-^11. 



Now, existence of a uniform Tc{N, a, K) follows from compactness of K by standard argu- 
ments. □ 

We also need a lemma which controls from below the growth of the norm of free 
solutions corresponding to rough eigenf unctions. We denote by Pp the spectral projection on 
the pure point spectrum of the operator L. 

Lemma 3.3. Assume that not a single eigenvector of the operator L belongs to H^{r). Let 
K C S be a compact set. Consider the set Ki = {(f) E K \ ||Pp0|| > 1/2}. Then for any B > 
we can find Np{B,K) and Tp{B,K) such that for any N > Np{B,K), any T > Tp{B,K) 
and any (j) E Ki, we have 

T 

i j \\PNe^'^'Pp<P\\ldt>B. (3.2) 



Remark. Note that unlike in (j3.1|) . we have the norm in (j3.2j) . 

Proof. The set Ki may be empty, in which case there is nothing to prove. Otherwise, let us 
denote by Ej the eigenvalues of L (distinct, without repetitions) and by Qj the orthogonal 
projection on the space spanned by the eigenfunctions corresponding to Ej. First, let us show 
that for any P > there is N{B, K) such that for any (p ^ Ki we have 

J2\\PNQjm>2B (3.3) 

j 

if N > N{B, K). It is clear that for each fixed (p with PpCp 7^ we can find N{B, cp) so that 
()3.3|) holds, since by assumption QjCp does not belong to whenever Qj(p 7^ 0. Assume that 
N{B,K) cannot be chosen uniformly for G i^i. This means that for any n, there exists 
(pn G Ki such that 

J2\\PnQ,(Pn\\l<2B. 
3 
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Since Ki is compact, we can find a subsequence ni such that converges to G -ft'i in if 
as — > oo. For any and any n/^ > N we have 



The last inequahty follows by Fatou's Lemma from the convergence of (pm to (f) in H and the 

1 /2 

fact that llPHj^Qj-^Alli < A„;^ for any n;^. But now the expression on the right hand side 
is less than or equal to 

liminfV||P„,Q,0„J|2<25. 



Thus ^ ||P/vQj(/'||i < 25 for any A^, a contradiction since G -ft'i. 
Next, take (p E Ki and consider 



^ / \\PNe''^%mdt = Y, 



^i{Ej-Ei)T _ ^ 



{rP^Q,<j),P^Qi<P). 



(3.4) 



In ()3.4|) . we set (e**^^^ ^'•''^ ~ 1)/ ^(-E'j ^ -f'/)^ = 1 if j = /. Notice that the sum above converges 
absolutely. Indeed, P/vQj^ = Yl!i=i{Qj^i^i)^ii ^^"^ (y^ii^k) = \^ik-, therefore 

AT 

{TPNQj4>,PNQi(p) = j4>, ej) {Qi(j), Cj) . 

i=l 

Therefore the sum on the right hand side of ()3.4|) does not exceed 

N N 

5^ A, 5^ I (g,-0, e,) (Q,0, e,) I < A;v 5^ 5^ ||Q^0|| I , e,,) (^^^^ 

i=l j,l 1=1 j,l 

N 

^ A^5^5^||g;</.|n(g,</./||g,0||,e,)|2 < A^iv, (3.5) 



with the second step obtained using the Cauchy-Schwartz inequality, and the third by 
||ei|| = 1. Then for each fixed N, we have by the dominated convergence theorem that the 
expression in (j3.4j) converges to J2j l|r"^''^-P/vgj0||^ = ||-P/vgj0|li as T — > cxd. Now assume 
N > Np{B,K) = N{B,K), so that (jSISl) holds. We claim that we can choose Tp{B,K) so 
that for any T > Tp{B, K) we have 



i{E,-Ei)T _ I 

Y^-r^ — ^(rp^g,0, p^Qi^ 



<B (3.6) 



for all G -f^i. Indeed, this follows from convergence to zero for each individual as T — > oo, 
compactness of Ki, and uniform continuity of the expression in the middle of ()3.6|) in for 
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each T (with constants independent of T) . The latter is proved by estimating the difference 
of these expressions for (pjtp E Ki and any T by 

Y^K^PnQA, PNQii<P -m + \(^PnQj{<P - V'), PNQii^)\, 

which is then bounded by 2XNN\\(f) — using the trick from ()3.5j) . Combining ()3.3|1 and 
()3.(ij) proves the lemma. □ 

We can now complete the proof of Theorem I1.41 

Proof of Theorem \1.4\ Recall that given any r, 5 > 0, we need to show the existence of eo > 
such that if e < eo, then ||0^(T/e)|| < S for any initial datum 0o ^ H, \\4>o\\ = 1- Here 
is the solution of ()2.4p . Let us outline the idea of the proof. Lemma 12.31 tells us that if 
the norm of the solution is large, relaxation is happening quickly. If, on the other 
hand, ||0^(to)||i < Aa/||0^(to)|P, where M is to be chosen depending on t,6, then the set of 
all unit vectors satisfying this inequality is compact, and so we can apply Lemma and 
Lemma 13.31 Using these lemmas, we will show that even if the norm is small at some 
moment of time tq, it will be large on the average in some time interval after tq. Enhanced 
relaxation will follow. 

We now provide the details. Since T is an unbounded positive operator with a discrete 
spectrum, we know that its eigenvalues A„ ^ cxo as n — > cxd. Let us choose M large enough, 
so that e-^*^^/^° < 6. Define the sets K = {(p e S \ < Xm} C S and as before, 
Ki = {(f) E K \ \\Pp(f)\\ > 1/2}. It is easy to see that K is compact. Choose N so that N > M 
and > Np{5XM, K) from Lemma f3. 31 Define 

n = max |Tp(5Am, K),T,{N, ^fc, ^) } , 

where Tp is from Lemma 13.31 and Tc from Lemma 13.21 Finally, choose eo > so that 
Ti < r/2eo, and 

n 


where B{t) is the function from condition ()2.7|1 . 

Take any e < eo. If we have ||0^(s) ||^ > Am for all s G [0, r] then Lemma ITSl implies 
that ||0^(T/e)|| < e"^^""^ < ^ by the choice of M and we are done. Otherwise, let tq be the 
first time in the interval [0, r/e] such that ||0'^(to)||^ < Am||0''(to)|P (it may be that To = 0, of 
course). We claim that the following estimate holds for the decay of 110*^(^)11 on the interval 
[t-o,t-o + Ti]: 

||0^(ro + n)f <e-^--^/2°||0^(ro)f. (3.8) 
For the sake of transparency, henceforth we will denote 4''^{tq) = (pQ. On the interval 
[ro,ro + Ti], consider the function (f)^{t) satisfying -j^cjPit) = iLcjj^it)^ <P^{tq) = 4>q- Note that 
by the choice of eo, (j3.7p and Corollarv 12.51 we have 

||0^(t)-/(t)|p<-^||0oir (3.9) 
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for all t G [to,tq + ti]. Split 0°(t) = </>c(^) + 'Ppi't)^ where 0c,p also solve the free equation 
^(pc,p(t) = iL(f)c,p{t), but with initial data Pc0o and Pp4>o at t = tq, respectively. We will now 
consider two cases. 

Case I. Assume that ||Pc0o||^ > fll^oP; or, equivalently, ||Pp0o|P ^ ill0oP- Note that 
since 0o/||0o|| & K hy the hypothesis, we can apply Lemma Our choice of ri implies 
that 

To+n 

- f WPNMtWdtK^Uof- (3.10) 

TO 

By elementary considerations, 

ii(/-p^)0°(t)f > l\\{i-PN)Mt)r-\\{i-PN)Mt)r > lmt)r-l\\PNMt)r-\\Mt)r- 

Taking into account the fact that the free evolution e*^* is unitary, Aat > Am, our assumptions 
on ||Pc,p0o|| and ()3.10p . we obtain 

TO+ri 

^ I " ^ /.0/-MI2 J+ ^ 1 II J. ||2 



„J-P^)/(t)||^dt> (3.11) 

Ti y 10 



TO 



Using ()3.9p . we conclude that 

TQ+Tl 

i 1 \\{i-p^mt)rdt>^\\4>or- (3.12) 



TO 

This estimate implies that 

TO+Tl 



ldt>^\\M'- (3.13) 

Combining ()3.13j) with ()2.6|) yields 

Il0^(ro + ri)f < (1-^) ||0^(ro)f <e-W^°||0^(ro)ir. (3.14) 

This finishes the proof of ()3.8|1 in the first case since Aat > Am- 

Case II. Now suppose that ||Pp0o|P > ill0o|p. In this case (/)o/||</'o|| G Ki, and we can 
apply Lemma f3. 31 In particular, by the choice of and ri, we have 

TO+Tl 

- / \\PN(t>pmldt>5XMUof. (3.15) 

Tl J 

TO 



DIFFUSION AND MIXING IN FLUID FLOW 13 

Since f|3.1Up still holds because of our choice of Tq and ri, it follows that 

TO+Tl 

I \\PNMmidt<^Uor- (3.16) 



TO 



Note that the if-norm in ()3.10jl has been replaced in ()3.16j) by the iJ^-norm at the expense 
of the factor of Ajv- Together, (|3.15p and (j3.16p imply 



TO+ri 



;i J \\PN<l>\t)\\ldt>2XMHof. (3.17) 

Finally, (HTTTI) and dHS)) give 



TO 



\PNrmidt>^Hor (3.18) 



TO 

since ||P/v0' - Pn<I)°\\1 < Ajv||0' - • As before, (j3T8j) implies 

||0^(ro + ri)ir<e-^--Ml0^(ro)ir, (3.19) 

finishing the proof of ()3.8|) in the second case. 

Summarizing, we see that if ||0^(to)||i < AA/||0^(ro)|p, then 

||0^(ro + ri)f < e-^-^-^/2°||0^(ro)f . (3.20) 

On the other hand, for any interval / = [a, b] such that ||0'^(t)||i > Aa/ 110*^(^)11^ on J, we have 
by Lemma f2. 31 that 

Wib)f < e-2^'>^^('-")||0^(a)f . (3.21) 

Combining all the decay factors gained from ()3.20j) and ()3.21|) . and using Ti < r/2e, we find 
that there is T2 G [r/2e, r/e] such that 

||0'(T2)ir < e-^^^^"2/2° < e-^'^^"/^" < 6^ 

by our choice of M. Then fl2.6|) gives ||0'^(T/e)|| < ||0'^(t2)|| < 6, finishing the proof of 
Theorem 11.41 □ 

4. Examples With Rough Eigenvectors 

It is not immediately obvious that condition ()2.7p . || 6*^*0 ||i < i?(t)||0||i for any 00 G H^, 
is consistent with the existence of eigenvectors of L which are not in H^. The purpose of 
this section is to show that, in general, rough eigenvectors may indeed be present under the 
conditions of Theorem 11.41 We provide here two simple examples of operators F and L in 
which ()2.7p is satisfied and L has only rough eigenf unctions. In both cases L will be a discrete 
Schrodinger operator on Z+ resp., more generally, a Jacobi matrix, and F a multiplication 
operator. One more example with rough eigenfunctions will deal with an actual fluid flow 
and will be discussed in Section IHl 
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The first is an explicit example with one rough eigenvector that is a discrete version of the 
celebrated Wigner-von Neumann construction of an imbedded eigenvalue of a Schrodinger 
operator with a decaying potential. The second example is implicit, its existence being 
guaranteed by a result of Killip and Simon , and demonstrates that all eigenvectors of L 
can be rough while at the same time the eigenvalues can be dense in the spectrum of L. 

Example 1. Let F be the operator of multiplication by n on /^(Z+), Z+ = {1,2, . . . , }. 
Furthermore, let L be the discrete Schrodinger operator on /^(Z+) 

LUn = Un+l + Un-1 + VnUn 

for n > 1, with the potential 

■^-^ n even, 

n+2 ' 



2 

I—; 

-1 n = h 



= < ^IZT ^ > 1 odd. 



and the self-adjoint boundary condition mq = 0. Then L has eigenvalue zero with eigenfunc- 
tion u given by 

_ _ (-1)" 

U2n-l — U2n — — - — 

for n > 1, because then Lu = and u G £^(Z^). Note that u does not belong to if^(F). 

It is not difficult to show that L has no more eigenvalues in its essential spectrum [—2, 2] 
(for example, using the so-called EFGP transform, see for more details). The eigenvalue 
zero is a consequence of a resonant structure of the potential which is tuned to this energy. 
There may be (and there are) other eigenvalues outside [—2,2], with eigenfunctions that 
are exponentially decaying and so do belong to H^{r). It is also known that L has no 
singular continuous spectrum and it has absolutely continuous spectrum that fills [—2,2]. 
More precisely, the absolutely continuous part of the spectral measure gives positive weight 
to any set of positive Lebesgue measure lying in [—2,2] (see, e.g., |2S])- 

To get an example where we have only rough eigenfunctions, we will project away the 
eigenfunctions lying in H^. Namely, denote by D the subspace spanned by all eigenfunctions 
of L, with the exception of u. Denote P the projection on the orthogonal complement of D, 
and set F = PTP, L = PLP. Then F, L are self-adjoint on the infinite dimensional Hilbert 
space H = Pl^{Z~^), and by construction L has absolutely continuous spectrum filling [—2, 2] 
as well as a single eigenvalue equal to zero. The corresponding eigenfunction is u and it does 
not belong to H^(T) because 

(Fm, u) = {PTPu, u) = {Tu, u)>Y^ \n\ {n-^f = oo. 

n 

Let us check the conditions of Theorem 11.41 First, F is positive because F is. It is also 
unbounded and has a discrete spectrum. Indeed, let Hr C H he the subspace of all vectors 
(f) E H such that 

(F0,0)<i?(0,0). (4.1) 
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Then for each such we also have ()4.H) with F instead of F. By the minimax principle for 
self-adjoint operators this implies that A„ > = n, where A^, A^ are the n-th eigenvalues 
of F and F, respectively (counting multiplicities). 

Also, L is a bounded operator on H (since L is) and so ()2.H) is satisfied automatically. 
Finally, observe that for any G H^{r), we have 

|(FL0,0)| = \{PfPLP<P,<P)\ = |(fL0,0)| < CUWl (4.2) 

The second equality in ()4.2|) follows from the fact that L and P commute by construction and 
P(f) = for G The inequality in ()4.2|) holds since ||-f/0||i < C||0||i, which follows from 
the fact that L is tridiagonal and both A„+i/A„ and f„ are bounded. Now given G H^(r), 
set 0(t) = e*"0. Then 

<2|(FL0(t),0(t))|<C||0(t)||? 

by ()4.2|) . This a priori estimate and Gronwall's Lemma allow one to conclude that ()2.7p holds 
with B{t) = e*"*/^. This concludes our first example. 

Example 2. We let H = £^(Z"^) and define F to be the multiplication by e". In order to 
provide an example with a much richer set of rough eigenfunctions, we will now consider L 
to be a Jacobi matrix 

with a„ > 0, f„ G M and boundary condition uq = 0. We choose to be a pure point 
measure of total mass ^, whose mass points are contained and dense in (—2, 2), and define the 
probability measure dfj,{x) = dvi^x) + \x[-2,2]{x)dx. By the Killip-Simon ^Sj characterization 
of spectral measures of Jacobi matrices that are Hilbert-Schmidt perturbations of the free 
half-line Schrodinger operator (with a„ = 1, f„ = 0), there is a unique Jacobi matrix L such 
that a„ — 1, f„ G £^(Z+) and its spectral measure is /x. In particular, the eigenvalues of L are 
dense in its spectrum a{L) = [—2,2]. 

The conditions of Theorem 11.41 are again satisfied, with the key estimate 11-^^0 ||i < C'||0||i 
holding because A„+i/A„, a„, t>„ are bounded. Moreover, it is easy to show (see below) that 
the fact that eigenvalues of L are inside (—2, 2) and a„ — 1, Vn G i"^ imply that eigenfunctions 
of L decay slower than e"*"^ for some C. More precisely, if u is an eigenfunction of L, then 
limji(ii^ -|- M^_i)e*"^ = oo, and so obviously u ^ H^(r) (actually, u ^ H^{T) for any s > 0). 

To obtain the well-known bound on the eigenfunction decay, let u be an eigenfunction of 
L corresponding to eigenvalue E G (—2,2), that is, 

EUn = anUn+l + a^-l^n-l + ^^n^n (4.3) 

for n>l. Define the square of the Priifer amplitude of u by 

Rn = ul + - EUnUn-l = ^^^^ (w^ + M^i) + ^("n " ^n-l)^ > 
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and c„ = |a„ — 1| + |a„_i — 1| + \vn\ G After expressing Un+i in terms of Un and Un-i using 
(14 .31) . one obtains (with each |0(c„)| < CeCu) 



Rn+l , ^^ ^((1 + 0{Cn))ul + (1 + 0{Cn))ul_^) + ^(1 + 0(c„)) K - 

2 



;i+o(c„))- 



if E 7^ and 

(1 + 0{cn))ul + (1 + 0(c„))nLi + 0( 
if ii^ = 0. In either case, Rn+i/Rn = 1 + 0{cn), which means that 

n n 

Rn > Rno n " ^ECk) > Rno exp{-2CE ^ Cfe) > Rno exp(-2C£; || Cfe || 2 ^/n) 
fc=no+l k=no+l 

if riQ is chosen so that CECk < for k > uq. But then the definition of Rn shows that 
Umji(u^ + M^_i)e*"^ = oo for some C < oo. This concludes the example. 
We have thus proved 

Theorem 4.1. There exist a self- adjoint, positive, unbounded operator T with a discrete 
spectrum and a self-adjoint operator L such that the following conditions are satisfied. 

• IIL0II < and ||e^^*0||i < B{t)\\(p\\i for some C < oo, B{t) e L^^JO, oo) and any 

• L has eigenvectors but not a single one belongs to H'^{r) for any s > 0. 

Later we will discuss examples of relaxation enhancing flows on manifolds. One of our 
examples is derived from a construction going back to Kolmogorov jSHj , and yields a smooth 
flow with discrete spectrum and rough eigenfunctions. This example is even more striking 
than the ones we discussed here since the spectrum is discrete. However, the construction is 
more technical and is postponed till Section IHl 



5. The Fluid Flow Theorem 

In this section we discuss applications of the general criterion to various situations involving 
diffusion in a fluid flow. First, we are going to prove Theorem 11.21 Most of the results we 
need regarding the evolution generated by incompressible flows are well-known and can be 
found, for example, in jSU] in the Euclidean space case. There are no essential changes in the 
more general manifold setting. 

Proof of Theorem ] 1.^ It is well known that the Laplace-Beltrami operator A on a compact 
smooth Riemannian manifold is self-adjoint, non-positive, unbounded, and has a discrete 
spectrum (see e.g. |^). Moreover, it is negative when considered on the invariant subspace 
of mean zero functions. Henceforth, this will be our Hilbert space: H = L'^{M) 1. 
Obviously it is sufficient to prove Theorem 11.21 for (po & H (i.e., when = 0). The Lipschitz 
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class divergence free vector field u generates a volume measure preserving transformation 
defined by 

= u(<l>t(x)), $o(x)=x (5.1) 

(see, e.g. [30] )■ The existence and uniqueness of solutions to the system ()5.1|) follows from 
the well-known theorems on existence and uniqueness of solutions to first order systems 
of ODEs involving Lipschitz class functions. With this transformation we can associate a 
unitary evolution group in L'^{M) where U^f{x) = /($_t(a;)). It is easy to see that H is an 
invariant subspace for this group. The group f/* corresponds to e^^* in the abstract setting of 
Section H Since |(t/*/) = -u ■ V(f/*/) for all / G H\M) (the usual Sobolev space on M), 
we see that the group's self-adjoint generator, L, is defined by L = m ■ V on functions from 
H^{M). It is clear that condition (EH) is satisfied, since ||m ■ V/|| < C||/||i for all / G if^ 
It remains to check that the condition ()2.7|) is satisfied, that is, ||e*^*/||i < -B(i(:)||/||i. Notice 
that if u{x) is Lipschitz, so is $t(x) for any t. This follows from the estimate (in the local 
coordinates and for a sufficiently small time t) 



|$t(x) - ^t{.y)\ <\x-y\ + J |u($,(x)) - u{^siy)\ ds. 



Applying Gronwall's lemma, we get 

|$t(x) - ^tiy)\ <\x- ?/|ell''l'L'p* 

for any x, y. Now by the well-known results on change of variables in Sobolev functions (see 
e.g. inZ!) and by the fact that $f is measure preserving, we have that 

||f/Vl|l<Cl|$,||Lipl|/l|l. 

This is exactly ()2.7j) . and the application of Theorem II .41 finishes the proof. □ 

The criterion of Theorem 11.41 can be applied to boundary value problems as well. For 
the sake of simplicity, consider a bounded domain f2 C M'^ with a boundary dQ. Let 
u G Lip(f2) be a Lipschitz incompressible flow such that its normal component is zero on the 
boundary: u{x) ■ h{x) = for x G dQ, with h{x) the outer normal at x. Let (f)^{x,t) be the 
solution of 

dt(f)^(x, t) + Au- V0^(x, t) - A^^(x, t) = 0, 0^(x, 0) = 0o(x), = if x e dQ, (5.2) 

on 

where the Neumann boundary condition is satisfled in the trace sense for almost every t > 0. 
The existence of solution to ()5.2|) can be proved similarly to Theorem 12.11 

Theorem 5.1. In the Neumann boundary conditions setting, the flow u G Lip(i7) is re- 
laxation enhancing according to the Definition if and only if the operator u ■ V has no 
eigenfunctions in H^{Q) other than the constant function. 

Proof. The proof is essentially identical to that of Theorem 11.21 The Laplacian operator 
with Neumann boundary conditions restricted to mean zero functions plays a role of the self- 
adjoint operator F. The condition u-n = ensures that the vector fleld u generates a measure 
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preserving flow via (jS-H) . and thus the corresponding evolution group is unitary. The 

estimates necessary for Theorem 11.41 to apply are verified in the same way as in the proof of 
Theorem 11.21 □ 

To treat other types of boundary conditions, such as Dirichlet, one needs to modify the 
relaxation enhancement definition. This is due to the fact that in this case the solution of 
fjl.lj) always tends to zero, rather than to the average of the initial datum. 

Definition 5.2. Let (j)^{x, t) solve evolution equation ()5.2j) . hut with Dirichlet or more general 
heat loss type boundary conditions 

-^{x,t)+aix)<i)^ix,t) = Q, xedQ, a(x)eC(dQ), (t(x) > (5.3) 
on 

where n is the outer normal to dQ. Then we call the divergence free flow u G Lip(i7) relaxation 
enhancing if for every t and 5 there exists A{t,5) such that for A > A{t,6) and ||0o||L2(n) = 1 
we have ||0^(a;, r)||i2(n) < 6. 

Remarks. 1. Note that cr(a;) = oo is not excluded and may lead to the Dirichlet boundary 
conditions on a part of the boundary. 

2. The more general definition encompassing both Definitions 11.11 and 15. 21 would assume that 
the solution tends to a certain limit and define relaxation enhancement in terms of speed up 
in reaching this limit. 

It is well known that the Laplace operator with boundary conditions (j5.3p is self-adjoint 
on the domain of H'^{Q) functions satisfying (j5.3p in the trace sense in L'^{dfl). We denote 
this operator A^-. The corresponding Hl{Q) space is the domain of the quadratic form of A^-, 
consisting of all functions e H^{fl) such that Jg^a{x)\(j){x)\'^ds is finite. In the Dirichlet 
boundary condition case, formally corresponding to a{x) = oo, we obtain the standard space 
Hq[Q). Then we have 

Theorem 5.3. In the case of the heat loss boundary condition (j5.3|) . the flow u G Lip(f2) 
satisfying u ■ h = on the boundary is relaxation enhancing according to Definition \5.^ if and 
only if the operator u ■ V has no eigenfunctions in H^{Q). 

Proof. In the case of heat loss boundary conditions, it is well-known that the principal eigen- 
value of Ao- is positive, so we can set T = —A„. Our space H is now equal to L'^{Q). The rest 
of the proof remains the same as in Theorem 15.11 □ 

We note that the case of Dirichlet boundary conditions has been treated in j2j in a more 
general setting u G L'^{Q) and without the assumption u ■ fi = 0. The methods of [2] 
are completely different from ours, and rely on the estimates on the principal eigenvalue of 
—A + Au ■ V and positivity of the corresponding eigenfunction. In particular, as described 
in the introduction, these methods do not seem to be directly applicable to the study of 
the enhanced relaxation in the case of a compact manifold without boundary or Neumann 
boundary conditions, where the principal eigenvalue is always zero. The results of j2] show 
that in the Dirichlet boundary condition case, the flow u is relaxation enhancing in the sense 
of Definition 15.21 if and only if u does not have a first integral in Hq{Q). In other words, if 
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and only if the operator u ■ V does not have an Hq{Q) eigenfunction corresponding to the 
eigenvalue zero. The discrepancy between this result and Theorem 15.31 mav seem surprising, 
but in fact the explanation is simple. 

Proposition 5.4. Let u G Lip(f2). If 4> & H^{Q) (H^{Q)) is an eigenfunction of the operator 
u ■ V corresponding to the eigenvalue iX, then |0| G H^{Q) (H^{Q) ) and it is the first integral 
of u, that is, u ■ V|(/>| = 0. 

Proof. The fact that \(f)\ G follows from the well-known properties of Sobolev functions 
(see e.g. [8 ). A direct computation using u ■ V0 = iX(j) then verifies that u ■ =0. □ 

As a consequence, when a ^ 0, the condition of no eigenfunctions in the statement of 
Theorem 15.31 can be replaced by the condition of no first integrals in H^. In the settings of 
Theorems 11.21 and 15. ![ the above argument still applies but does not allow to change their 
statements. Indeed — on one hand the operator u ■ V always has eigenvalue zero with an 
eigenfunction that is smooth, namely a constant. Existence of this first integral, however, 
tells us nothing about relaxation enhancement. On the other hand, existence of mean zero 
eigenfunctions need not guarantee the existence of a mean zero first integral, as can be 
seen in the following well-known example. 

Example. Let M = T'' be the flat d-dimensional torus with period one. Let a he a. d 
dimensional constant vector generating irrational rotation on the torus (that is, we assume 
that components of a are independent over the field of rationals). It is well known that 
the flow generated by the constant vector field a is ergodic but not weakly mixing. The 
self-adjoint operator L = ia -V has eigenvalues 27rQ; ■ k, where k are all possible vectors with 
integer components. The corresponding eigenfunctions are e"^'^*'^'^, x G T'^. Their absolute 
value is 1, which is a first integral of a, but there are no other first integrals. In particular, 
every non-constant eigenfunction of L corresponds to a non-zero eigenvalue. Thus, this flow is 
not relaxation enhancing even though it has no first integrals other than a constant function. 



Finally, we show that the norm in the Definitions II. ![ 15.21 can be replaced by other 
norms with 1 < p < oo without any change to the statements of Theorems 11.21 I^TTl 15.31 This 
result is important for applications to quenching in reaction-diffusion equations. 

Theorem 5.5. Theorems M.iA f3~?l 15.^ remain true if, in Definitions M.IV \5.'A "H0o||l2 = 
1" IS replaced by "||0oIUp = 1" and 1|0^(r) - 0||l2 < 6" (resp. "||0^(r)||L2 < 5") by 
1|0^(r)-0|U, <5" (resp. "\\<P^{t)\\l. < S") for any p,q e [l,oo]. 

For the sake of consistency of notation, we will consider the compact manifold case. The 
case of a domain Q with Dirichlet or heat-loss boundary conditions is handled similarly (see 
below) . 

We start with the proof of a general L°° estimate for solutions of 

ipt + v - Vij - Atp = (5.4) 

on a compact manifold M. The point is that this estimate will be independent of the 
incompressible flow v and so, in particular, of the amplitude A in ()1.1|) . It appeared, for 
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= -2||V0||i. < -2CM?v^'UU < -2C,U&Ull (5.7) 



example, in J^Ij where the domain was a strip in M?. The crucial ingredient of the proof 
was a Nash inequality. In the general case, we follow a part of the argument, but our proof 
of the corresponding inequality ()5.fij) is different. 

Lemma 5.6. For any smooth Riemannian manifold M of dimension d and any £ > (resp. 
e > 0) if d > 3 (resp. d = 2), there exists C = C{M,e) > such that for any incompressible 
flow V G Lip(M) and any mean zero 0o ^ L'^{M), the solution of ()5.4|1 and 0) = (poix) 
satisfies 

mX,t)\\L^{M) < Cr"/'-^||0o||Li{M). (5.5) 

Proof. First note that by Holder and Poincare inequalities we have for any mean zero ijj G 
H^{M), any p > (c/ + 2)/4 (if d = 2, then for any p > {d + 2)/4), and some Cp, 

That is, 

l|V^||i.>C,||^||J^t^||^||-? (5.6) 

for q = 2/{2p - 1) so that g < 4/(i if d > 3 and g < 2 if d = 2. 

After multiplying ()5.4j) by (p and integrating over M we obtain for t > 

d 
dt 

The last inequality follows from the positivity of and the preservation of norms of solutions 
of ()5.4|) with initial conditions 0o,± = max{±0o, 0}, which shows that H'/'IIli is non-increasing. 

Next we divide ()5.7p by — ||0||^t^ and integrate in time to obtain ||0(x, t)||^| > qCqt\\(j)o\\~i . 
This in turn gives with a new Cg 

mx,t)\\L^<Cgt-'/^4>o\\L^. (5.8) 

Hence we have showed that ||'P((f) || 2,1^^2 < Cqt~^/'^ where Vt{v) is the solution operator for 
()5.4j) . But since Vtiy) = {Vt{—v))* is the adjoint of the operator Vt{—v), which satisfies the 
same bound, we obtain 

\\V2t{v)\\Li^L^ < \\Vt{v)\\Li^L4'Pt{v)\\L^-,Lo^ = \\Vtiv)\\Li^L4'Pt{-v)\\L^^L^ < C'gt-^^" 

which is ()5.5|) . □ 

Proof of Theorem \5.^ Assume for simplicity that the total volume of M is equal to one. 
Then it is clear that Lemma [5.61 also holds with 

||0(x,t)l|LP(M) < Ct-'^/'-'UohHM) (5.9) 

in place of (|5.5|) . with any p,q E [1, 00] and the same C. 

Assume now that we know that for some u and p,q > 1, given any t,S > 0, we can find 
Apg{T,6) such that ||0^(x, r)||LP < S for any A > Ap^g{T,6) and any mean zero 0o ^ L'^{M) 
with 1 1 00 1 = 1. Take any other p',q' > 1. Then for any A > Ap,g(T, S), 

||0^(x,3r)||^,, < Cr-^/2-^||0^(x,2r)|U. < 6Ct-'/'~4<P^{x,t)U. < 6{CT~'/''n'\\ML.' 
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This shows that when Ap q{T, 6) exists for some p, q and all r, 6, for any p', q', r, 6 we have 
Api^q'^T, 6) = Ap^qi^r/S, 6C^'^{t/3Y~^^^) and so Ap'^q'^r, 6) exists for all r, 6. That is, Definition 
11.11 describes the same class of flows regardless of which — > L"? decay it addresses. This 
finishes the proof. □ 

We note that in the case of a bounded domain Q with Dirichlet boundary conditions, the 
proof is identical. When we have heat-loss boundary conditions, the only change is that the 
equality in ()5.7|) reads 

|||0||i. = -2(||V</>||i. + ||aV20||i.(,^)) 
and the Poincare inequality is replaced by 

< CpiML^ + II VV^IUO < {Cp + Ao ^/')(|| V^IU^ + IW'/'^Wman)) 



which is due to the Sobolev inequality and the fact that the principal eigenvalue Aq of the 
Laplacian on Q with heat-loss boundary conditions is positive. 

6. Examples of relaxation enhancing flows 

Here we discuss examples of flows that are relaxation enhancing. Most of the results in this 
section are not new and are provided for illustration purposes. According to Theorem 11.21 a 
flow u G Lip(M) is relaxation enhancing if all of its eigenfunctions are not in H^{M). One 
natural class satisfying this condition is weakly mixing flows - for which the spectrum is 
purely continuous. Examples of weakly mixing flows on go back to von Neumann jSS] 
and Kolmogorov [2H] • The flow in von Neumann's example is continuous; in the construction 
suggested by Kolmogorov the flow is smooth. The technical details of the construction have 
been carried out in see also [201 ^ review. Recently, Fayad ^3] generalized this 
example to show that weakly mixing flows are generic in a certain sense. For more results on 
weakly mixing flows, see for example ^1 120] • To describe the result in ^2] in more detail, 
let us recall that a vector a in M.'^ is called /?-Diophatine if there exists a constant C such 
that for each G Z'^ \ {0} we have 

C 

inf I (a, k) + p\ > 



The vector a is Liouvillean if it is not Diophantine for any (3 > 0. The Liouvillean numbers 
(and vectors) are the ones which can be very well approximated by rationals. 

Example 1. Consider the flow on a torus T'^"'"^ that is a time change of a linear translation 
flow: 

with a smooth positive function F{x,y). Such flows have a unique invariant measure dfi = 
F{x,y)dxdy. Let us denote by C"'(T'^,M'^) the set of functions on the torus that are 
positive. We have 
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Proposition 6.1 ([13j). Assume the irrational vector a G M'^ is not (3-Diophantine, for some 
(3 G U +00. Then, for a dense Gs of functions F in C'^+°'(T°', M"*") the flow is weakly 
mixing (for the unique invariant measure F{x,y)dxdy). 

To obtain examples of relaxation enhancing flows, we can now consider the generic flows 
of Proposition 16. II on the torus with a metric such that the volume element is F{x,y)dxdy. 
Alternatively, we can just view a problem in a weighted space and consider the operator 
|;V(FV), which is self-adjoint on L'^{T'''~^^ , Fdxdy), instead of A. It is also straightforward 
to obtain more physical examples of the relaxation enhancing flows that are incompressible 
with respect to the usual flat metric. Indeed, assume for the sake of simplicity that we are 
working with a unit torus and that the total integral of F is also one. Then it is not difficult to 
construct an explicit measure preserving invertible transformation Z from T'^~^^{F{x, y)dxdy) 
to T'^'^^ (dpdq) , as smooth as the function F (for general results on existence of such maps see 
[SIl). If we denote w{x,y) = {a/ F{x,y), 1/F{x,y)) the vector field in ()6.1|) . then the vector 
field u{p,q) = Z o w o Z~^ is going to be incompressible. Moreover, the unitary evolutions 
generated by w and u in L'^{F{x,y)dxdy) and L'^{dpdq) respectively are unitary equivalent 
and so have the same spectra. 

We now describe an example of a different class of flows to which Theorem 11.21 applies. 
Namely, we will sketch a construction of a smooth incompressible flow u{p, g), V ■ m = 0, on 
a torus such that it has a purely discrete spectrum but none of the eigenfunctions are in 
iif^(T^). We could not flnd an exact statement regarding the existence of such flows in the 
literature, although the idea of the construction appears in j2H| and the result is presumably 
well known in the dynamical systems community. In particular, it follows in a fairly direct 
way for example from considerations in PQI^HI- We briefly sketch the construction, without 
presenting well known technical details. 

Example 2. Let us denote by $" the flow on the torus generated by u and by t/* the flow 
on L2(T2) generated by •I'J': (t/7)(x) = /(^'^^(x)). 

Proposition 6.2. There exists a smooth incompressible (with respect to the Lebesgue mea- 
sure) flow u{x,y) on a two-dimensional torus so that the corresponding unitary evolution 
f/* has a discrete spectrum on L^(T^) but none of the eigenfunctions of U are in H^lT"^). 

Proof. The example will be given by a flow of type ()6.1|) . with d = 2 and appropriately 
chosen a and F{x,y). We remark that while the form ()6.1|) may seem quite special, in fact 
any analytic flow in two dimensions with an integral invariant can be mapped analytically to 
the linear translation flow with some a, F (see e.g. |2H])- The idea of the construction 
is to flnd a smooth flow which can be mapped to a constant flow {a, 1) by a measure 
preserving map S with very low regularity properties. Since the eigenfunctions of the constant 
flow are explicitly computable, we can compute the eigenfunctions of the original flow. Due 
to the roughness of 5", these will prove highly irregular. To obtain an incompressible flow, we 
will then proceed as in the flrst example. 

In order to flnd such a flow w, we start with a smooth periodic function Q G C°°(§^) and 
an irrational number a G M so that the homology equation 

i?(e + «)-i?(0 = Q(0-i, (6.2) 
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has a solution R{^) that is very rough. Note that for ()6.2|) to have a measurable solution the 
function Q{C,) should satisfy the normalization [Ij 

Jo 

The following Proposition is a particular case of Theorem 4.5 of [SDj. 

Proposition 6.3. Let a be a Liouvillean irrational number. There exists a C°°(S"'^) function 
Q{^) so that the homology equation \6. iij) has a unique (up to an additive constant) measurable 
solution R{S) : §^ ^ M such that for any A G = M\{0}, the function Rx{C) = e*^-^(«) is 
discontinuous everywhere. 

Note that without loss of generality we may assume that Q{C,) is positive - otherwise we 
choose M so that Q{^)+M > 1 and consider arescaled function Qm{0 = 
Then the function Rm{0 = R{0/i^ + 1) is the solution of ()6.2|1 with Qm on the right side 
and, of course, Rm{0 has the same properties as -R(0- 

Given a Liouvillean irrational number a and a function (5(0 that satisfies the conclusion 
of Proposition 16. HI we define a function F{x,y) on the torus as follows. Choose m > so 
that m < mmQ{s) and a smooth function iply) > such that 

/ ij{y)dy=l. (6.3) 

and, in addition, 

i'iy) = for < y < yo and yi < y < I with yQ close to zero and yi close to one. (6.4) 
The choice of m ensures that the function 

F{x,y) = m + ■ilj{y){Q{x — ay) — m), 0<x,y<l (6.5) 

is positive - then we extend F{x, y) periodically in both variables to the whole plane M^. 
The resulting function is smooth because of ()6.4|) and, in addition, it has total mass equal to 
one. The normalization ()6.3|) implies that the functions F and Q are related by 

Q(0 = / F{^ + az,z)dz. (6.6) 
Jo 

Now, the required transformation S : {x,y) ^ (X, Y) is defined by [2H1 ISH] 

X(x, y) =x + a{Y{x, y)-y), Y{x, y) = T{x - ay, y) + R{x - ay) (6.7) 

with the function R{x) that satisfies ()6.2p . and T{x,y) defined by 

ry 

T{x,y) = / F{x + az,z)dz. 
Jo 

Note that the transformation ()6.7p implies that x — ay = X — aY, and so it preserves the flow 
trajectories. The homology equation ()6.2|) together with the definition ()6.5p of the function 
F{x, y) imply that S is well-defined as a mapping — > T^. It is also straightforward to check 
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that S maps the flow w onto the uniform flow Wunif = !)• One can also verify that S is 
invertible with measurable inverse, and is measure preserving: 

[S*f]{x,y)F{x,y)dxdy = J f{S{x,y))F{x,y)dxdy = j f{X,Y)dXdY (6.8) 

for any function / G C(T^). Hence, S* may be extended as an operator F^^dxdy) — > L'^{dfi) 
with the preservation of the corresponding norms. We conclude that the unitary evolutions 

and U^nif generated by the flow w given by 1)6.11) and the uniform flow Wunif, respectively, 
are conjugated by means of a unitary transformation S* : L"^ {T"^ , dX dY) L^(T^,(i/i) 
- we have Ul^^J = [S*]^^UljS* . It follows that and U!^nif have the same spectrum: 

= 2'n'ina + 27r2/, l,n E 7^. It also follows that the eigenfunctions of the operator Uy^ may 
be written as 

1p^l{x y) = g2™X(a;,j/)+27ri«y(a;,j/) _ ^2mn{x-ay+aY{x,y))+2nilY{x,y) g-j 
^2iTin(x~ay) ^{2nina+2iTil)(T(x~ay,y)+R{x~ay)) (^(x y^g{'^'''"*"'^+^''^*0-^{^'~'^J/) 

with a smooth function C{x,y) G C°°{[0, 1]^) (note that the function ({x,y) is not periodic 
in {x,y)). In order to verify that ip^i are not in H^{T^) it suffices to check that the function 

Qxix,y) = e'^^'^^-'^y^ = Rxix - ay) 

is not in H^{[0, 1]^) for any real A 7^ 0. The function Rx{s) is defined in the Proposition 16.31 
and is everywhere discontinuous. If Q\{x, y) were in -ff^([0, 1]^), it would force R\[s) to be in 
if^(S^) and hence continuous but this function is discontinuous everywhere. Therefore, the 
eigenfunctions ip^i cannot be in H^iT"^). 

Finally, to obtain an incompressible flow, we introduce a smooth transformation Z : (x, y) — ; 
(p, q) be setting 

/■X _ in - 

p= F(s)ds, g = - / F(x,z)dz, where F(x) = / F(x,z)dz. 
Jo F{x) Jo Jo 

It is immediate to verify that Z maps the measure dfi onto the Lebesgue measure dpdq. 
Hence, the evolution group generated by the image u{p,q) of the flow w{x,y) will have the 
same discrete spectrum as Uw In addition, the eigenfunctions ipl^i of Uw are the images of 
the eigenfunctions iplt^i of u under Z*: ip'^i = Z*'iplt^i = ip'^i o Z. As the functions ijj'^i are not in 
if^(T^) and the map Z is smooth, it follows that all the eigenfunctions of the incompressible 
flow m(p, q) are not in H^iT"^). This finishes the proof of Proposition 16.21 □ 

7. Quenching in Reaction-Diffusion Equations 

In this section we describe the application of our results to questions of quenching in 
reaction-diffusion-advection equations. We will consider the problem 

T^{x,t) + Au-VT^{x,t)-AT^{x,t)=f{T^{x,t)), T^(x, 0) = To(x) (7.1) 

on a smooth compact Riemannian manifold M with To{x) G [0, 1]. Here T is the (normalized) 
temperature of a premixed flammable gas that is advected by the incompressible flow Au{x) G 
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Lip(M). The nonlinear right hand side term accounts for temperature increase due to burning 
and will be assumed to be of ignition type. That is, 

(i) /(O) = /(I) = and f(T) is Lipschitz continuous on [0, 1], 

(ii) 3^0 e (0, 1) such that /(T) = for T G [0, ^o] and /(T) > for T G (^o, !)• ^'^'^^ 

This shows, in particular, that T remains in [0,1]. The main question will be under what 
conditions on u one can always choose A large enough so that for some time r > we have 
||T^(a;, r)llLoo(^) < 9q, that is, quenching — extinction of flames — happens. Of course 
this question is meaningless for certain initial data Tq. Namely, if ||To||^i > ^ovol(M) or 
ll^ollii = 6'oVol(M) but To ^ Oq, then it is easy to show using 

|i|T^|U.= J f{T^{x,t))dx (7.3) 

that 1 1 T"^ 1 1^1 must be strictly increasing with the limit equal to the volume of M. This 
motivates the following definition. 

Definition 7.1. We say that u is strongly quenching if for any nonlinearity f as in ()7.2|) . and 

any solution of ()7.H) with initial datum Tq{x) G [0,1] with ||To||2,i(Af) < 6'ovol(M), there 
exists A(To, /) such that if A > A{Tq, f), then for some r > one has \\T"^(x, r) ||l°°(m) < ^o- 

Then we have 

Theorem 7.2. An incompressible flow u G Lip(M) is strongly quenching if and only if it is 
relaxation enhancing. 

Proof. Assume that the volume of M is one. First, u is strongly quenching when it is 
relaxation enhancing. Indeed, assume u is relaxation enhancing and let c be the Lipshitz 
constant for / so that /(T) < cT. If (p^ solves jHH) with 0o = Tq, then T^(x, t) < e^*0^(x, t) 
by the comparison principle. But this means 

\\T^{x,t)\\loo < e^^||0^(x,r)|Uoo < e^^(fo + ||0^(x,r) - 0|Uoo) 

which can be made as close to To = ||To||2,i(Af) < Oq as we wish by taking small enough 
T,6 > and A > Ai^oo{t, from the proof of Theorem 15.51 Since / was arbitrary, it follows 
that u is strongly quenching. 

Hence, we are left with proving that u being strongly quenching implies that u is relaxation 
enhancing. Assume this is not the case, that is, there exists an ignition nonlinearity /, a 
mean zero (po G L'^{M) and r, 5 G (0, 1) such that for a\\ A < oo and t < t, the solution (p"^ of 
(HUD satisfies \\(j)^{x,t)\\Li > 6. We can assume without loss of generality that ||0o||l°° < 1- 

Let To = 6^1 +700 where 7 = | min{6'o, 1— 6'o} and 9i G {0o—/3, 9o) with /? = | min{75, t5k, Oq} 
and K = min{/(T) | T G [^o + ^, |(2 + ^o)]} > 0. Note that To(x) G |[^o, 1 + 29o]. 

For each t < t let Bt = {x \ (p^^x, t) > |}. Since (p^{x, t) is mean zero with norm more 
than S, and t)| < 1, we must have \Bt\ > | for each t <t (recall that M has volume 

one). Thus 

T^{x,t)>9^ + ^4>9o + f 
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for X e Bt. If for some t < r there is a set C Bt with \B[\ > | and T^(x,t) > |(2 + 6^0) 
for X G Bl, then since ip^ = 61 + 70"^ < by the comparison principle (because tp^ satisfies 
(HH) with V'^(a;,0) = Tq{x)) and < inf{To(x)} < i)^ < ||To||l- < |(1 + 2^0) by the 
maximum principle, 

\\T\x, t)|Ui > |[|(2 + ^^o) - 1(1 + 29,)] + U\x, > f + > ^0. 

This is a contradiction because then u cannot be strongly quenching by the argument before 
Definition I7.ll 

Therefore for each t <t there must be a set B[' C Bt such that \B^\ > | and T^(x, t) G 
[^0 + f , |(2 + ^0)] for X G 5^ But then /(r^(x, t)) > k for x G 5^', and so (|ISD gives 

||T'^(x,r)|Ui>ei + rf >^^o 
and we have a contradiction again. Hence u has to be relaxation enhancing. □ 

Applications of our results on relaxation enhancement to quenching on infinite domains 
(where front propagation can occur as well) will be considered elsewhere. 
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